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Main characteristics of stationary anisotropic Poisson processes of cylin- 
ders (dilated fc-dimensional flats) in d-dimensional Euclidean space are stud- 
ied. Explicit formulae for the capacity functional, the covariance function, 
the contact distribution function, the volume fraction, and the intensity of 
the surface area measure are given which can be used directly in applications. 
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1 Introduction 

Porous fiber materials find vast applications in modern material technologies. Their 
use ranges from light polymer-based non- woven materials, see Helfen et al (2003), to 



fiber-reinforced textile and fuel cells as in Mukherjee and Wang 
percolation 



(|2006|). Their porosity. 
It is 



acoustic absorption and liquid permeability are of special interest, 
known that these properties depend to a great extent on the microscopic structure of 
fibers, in particular, on the orientation of a typical fiber. If all directions of fibers are 
equiprobable one speaks of isotropy. Many materials are made by pressing an isotropic 
collection of fibers together thus producing strongly anisotropic structures. As examples. 



pressed non-woven materials used as an acoustic trim in car production, see Schladitz 
etal|(p006l), paper making process as in ICorte and Kallmesl (11962]) , and Molchanov et al 



( 1993 ) , and gas diffusion layers of fuel cells, here Mathias et al ( 2003 ) , Manke et al ( 2007 ) 



can be mentioned; see Figure [T] 
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Figure 1: Microscopic structure of the gas diffusion layer as used in polymer electrolyte 
membrane fuel cells (courtesy of the Centre for Solar Energy and Hydrogen 
Research, Ulm) 



To quantify this dependence between the physical and the geometric structural prop- 
erties of porous materials, their intrinsic volumes (sometimes also called Minkowski 
functionals or quermassintegrals) are used. More formally, porous fiber materials are 
usually modeled as homogeneous random closed sets described in Matheron (19751 and 



Serra (1982). The mean volume and surface area of such sets in an observation window 



averaged by the volume of the window are examples of intensities of intrinsic volumes 
which are treated in detail in this paper. 

The intention of this paper is to give formulae for cylinder processes which can be used 
directly in applications, which is also demonstrated in the optimization example. Thus 
the focus is on stationary Poisson processes which are the most common in applications. 
A rather theoretical analysis can be found in the recent paper by jHoffmann (20091, 
where formulae for the curvature measures of a more general non-stationary model of 
Poisson cylinder processes can be found. In Weil (19871 the model for cylinders as used 
in this paper is introduced, and curvature measures for different kinds of (not necessarily 
Poisson) point processes are calculated. As opposed to that, in this paper formulae for 
the covariance function, the contact distribution function, and a different approach for 
the calculation of the specific surface area of Poisson cylinder processes are worked out 
which have straightforward applied value. 

As a model for fiber materials shown in Figure [T] we consider anisotropic stationary 
cylinder processes as homogeneous Poisson point processes in the space of cylinders. 
Isotropic models of this kind (named also processes of "thick" fibers, lamellae, membranes 



or Poisson slices) have been studied in detail, cf. Matheron (19751, Serra (19821, Davy 
(19781, Ohser and Miicklich (2000). See Schneider ( 1987[ ) for further references. In the 
present paper, we generalize some of their results to the anisotropic case. 

After giving some preliminaries on cylinder processes (Section [2]), we obtain formulae 



for the capacity functional, covariance function and contact distribution function in Sec- 
tion |3] In Section [4j we prove the formulae for the intensity of the surface area measure 
of anisotropic stationary Poisson processes of cylinders. Formulae for the intensities of 
other intrinsic volumes can be found in the recent paper by Hoffmann (20091. In the 
last section, we show how the volume fraction of an anisotropic Poisson process of cylin- 
ders can be maximized under certain constraints. In the solution, we use the formulae 
obtained in previous sections. 

Since the formulae obtained in Sections |3]-[4] are rather complex, examples in the most 
interesting dimensions 2 and 3 are given, which can be directly used in applications. 



2 Cylinder Processes 

Let G{k, d) be the Grassmann manifold of all non-oriented A;- dimensional linear sub- 
spaces of M*^, and © the cj-algebra of Borel subsets of G{k, d) in its usual topology. Let C 
(/C) be the set of all compact (compact convex) non-empty sets in M'^. Denote by TZ the 
convex ring, i.e., the family of all finite unions of non-empty compact convex sets. We 
provide these sets with the Hausdorff metric, and denote the resulting Borel cr-algebra 
of 7^ by ^R. 

Denote by z^fc(-) the fc-dimensional Lebesgue measure in M*-', and by n''{-) the k- 
dimensional Hausdorff measure. For any set S* C M'^ denote by S'^ the linear subspace 
of the vectors which are orthogonal to all elements of S. For ^ being a /c-dimensional 
flat (i.e. a /c-dimensional linear subspace) we denote the {d — A;)-dimensional Lebesgue 
measure in ^-^ by Let (oJk) be the volume (surface area) of a unit A;-dimensional 

ball, respectively. 

For a convex set K ^ K, and x G M"' let p{K, x) be the unique point in K which is the 
closest to x. Then there exist measures ^k{K, •) on B{W^), ior k = 0,...,d with 

d 

i/d({x eKe Brio) : p{K, x) e B}) = J2 r'^^'^Kd^MK, B), 

where Ki © K2 = {ki + k2\ki £ Ki, k2 G K2}, and Br{o) is the ball of radius r centered 
in the origin o. Furthermore we define $fc(0, B) =Q for all B G B{W^). These measures 
are called curvature measures. Since they are locally determined, they can be extended 
to functions with locally polyconvex sets as first argument in such a way that they 
remain additive. One should remark that these generalized curvatures measures are not 



necessarily positive, but signed measures. For a detailed introduction, see Schneider and 



VM] ( |2008| ). The intrinsic volumes of K can be defined as total curvature measures 



V^{K) = <^k{K, M"^) for fc = 0, . . . , d. 

Following the approach introduced in Weil 



sum of a flat ^ G G{k, d) and a set K (Z ^ 



( 1987 1 , we define a cylinder as the Minkowski 



with K £ IZ. Note that K is not limited 



to sets with an associated point in the origin. The flat ^ is also called the direction 
space of ^ ® i^T and K is called the cross section or base. For a cylinder Z = © we 
define the functions L{Z) = £^ and K{Z) = K. Furthermore, define Zk as the set of all 
cylinders which have a /c-dimensional direction space and base in IZ. For the volume of 
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L(Z) 

the cross-section of the cyhnder we introduce the notation A{Z) = u^_^^ {K{Z)). By 
S{K) we denote the surface area of a set K. In the case of K being the cross-section of 
a cyhnder K (B L we shah use this notation for the surface area of K in the space 

We cah a measure if on Zk locally finite if ip{{Z £ Zk\Z n K ^ 0}) < cxd for aU K £ C. 
Let M.{Zk) be the set of locahy finite counting measures on Z^ supphed with the usual 
cj-algebra 9Jt. A point process H on Zi. which is a measurable mapping from a probability 
space ($1, J^, P) into {Ai{Zk),Tl) is called a cylinder process. Its distribution is given by 
the probability measure P^: 97t — t- [0, 1], Ph(') = ^ ")• ^ cylinder process is called 
Poisson if is Poisson distributed with mean h.{B) for some locally finite measure 

A on Zk and all Borel sets B C Z^, and H(i?2), • • • are independent for 

all disjoint Borel sets Bi, B2, ■ ■ ■ , Bn C Zk, and all n > 2, see details in [Schneider and| 



Weil (20081. The measure A is called the intensity measure of H. The Poisson cylinder 
process is called simple, if it has no multiple points. This is the case if and only if A is 
diffuse. For the rest of this paper we assume that H is a simple Poisson cylinder process. 



In this case, the union [/=; = "^zgeZ is a random closed set, see (Schneider and Weil 



2008 p. 96), where we denote by Z G H the cylinders Z in the support set of H. 

The cylinder process H is called stationary if its distribution is invariant with respect 
to translations in M'^ and isotropic if it is invariant w.r.t. rotations about the origin. Let 
Z^ be the set of all cylinders K®^ with ^ G G{k, d),K C and for which the midpoint 
of the circumsphere of K lies in the origin. 



Following |Weil| ( |l987| , we define i: {x, Z) ^ x + Z foi x e M.'^ and Z G Z^. If E is 
stationary, then a number A > and a probability measure 9 on Z^ exist such that 



A{i{AxC)) = x[ ui[l\A)e{dZ) 



for all Borel sets A C W^,C C Z^. Then A is called the intensity and 6 the shape 
distribution of S. 



As shown in (Konig and Schmidt 1991 p. 61), 6 can be decomposed further. Anal- 
ogously to i define j : x ^) i->- © ^ for ^> G 7^ and G G{k,d). Then there exist 
a probability measure a on © (directional distribution of H) and a probability kernel 
/?: X G{k,d) — >■ [0,1] for which /3(-,^) is concentrated on subsets of such that for 
arbitrary i2 G D'l and G G the equation 

e{j{RxG))= [ (3{R,Oa{dO (1) 
Jg 

holds. 



3 Capacity functional and related characteristics 



In this section, we calculate the capacity functional (cf. (Stoyan et al, 1995, p. 195)) for 



the union set [/= of the stationary Poisson process H of cylinders with A;-dimensional 
direction space introduced as above. As a corollary, explicit formulae for the volume 
fraction, the covariance function, and the contact distribution function of follow 
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easily. It is worth mentioning that the resulting formula ^ for the capacity functional 
generalizes the formula in (Serra 1982 pp. 572-573), given for Poisson slices in M'^, and 



a model with this capacity functional has already been proposed in (Matheron 1975 
p. 148). 

3.1 Capacity functional 

For any random closed set X, the capacity functional Tx{B) = P(X f] B ^ fli), B ^ C, 
determines uniquely the distribution of X. 

Let TTjj{B) be the orthogonal projection of a set 2? C M*^ along a linear subspace ry C M"^ 
onto T]-^. 

Lemma 1. The capacity functional of the union set Us of the cylinder process H is given 
by 



Tu-SB) = 1 - exp |-A^^ v^^^:X\-K{Z) © ^z.(Z)(i?)) ^(dZ)| . 



(2) 



Proof. Let S be a compact set in M . Then by Fubini's theorem and (Schneider and 



Weil 2008 p. 96), we get 



1 - Tu^{B) = exp{-A({Z € Zk\Z DB^ 0})} = exp 1^- J 1{Z f] B ^ 0}A(dZ)| 

/ / l{(Z + x) n5 / 0}dx0(dZ) I 
/ ^ / 1{{K{Z) + x) n TTLiz) (B) + 0} dx ^(dZ) I , 



exp 

exp < — A 



where Z = x + Z . 

One can easily see that K{Z)+x hits t^l{z){B) if and only if x belongs to the Minkowski 
sum of —K[Z) and ttl{z){B). 
Thus we have 

1 - Tu^{B) = exp I -A J^^ j^^^^^ n{K{Z) + x) n ttl{Z){B) / 0} dx^(dZ) | 
= exp I -A j^^^^^ l{x G -K{Z) © TTLiz){B)}dxe{dZ)^ 
= exp |-A j^^ v^,[l\-K{Z) © TiLiz){B)) ^(dZ)| . 

□ 

A few remarks are in order. 
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It follows from the local finiteness of A that 



J^^ v^l\-K{Z) ® 7:Liz){B)) e{AZ) < oo, 



(3) 



cf. (Schneider and Weil 2008 p. 96, Theorem 3.6.3. and remark). 



• The choice of A; = yields the capacity functional of the stationary Boolean Model 
S' with the primary grain K and intensity A, cf. (Matheron 



1975 



p. 62): 



Ts'iB) = 1 - e 



-\Eui{-K(BB) 



Another important special case is that of K being a.s. a point. Then the model 



coincides with a /c-flat process H", cf. (Matheron 1975, p. 67) with the capacity 
functional 

T^.{B) = 1 - exp |-A j^^ y^^[l\^Liz){B)) 0(dZ)| . 

• The case B = {o} yields the volume fraction p = P(o G [/=) = 
-Ez^d(f/sn [0,1]^^) of [/=: 



p = Tu^{{o}) = 1 - exp I - A j^^ A{Z) e{dZ) | . 



(4) 



A generalization of this formula can also be found in Hoffmann ( 2009 1 in the non- 
stationary setting. 

Throughout this paper, we assume that f > 0, i.e. /^o A{Z)9{dZ) > 0. Thus, we 

. . k 

have p G (0, 1), cf. inequality ([3]). 
3.2 Covariance function 

In the following we investigate the covariance function of It is defined as 



Cus(h) = P(o,h e H), h G M'^, cf. (Stoyan et al 



1995 



p. 68). 



Because of the relation Cus{h) = P(o, h € H) = 2p—Tu^{{o, h}) it is closely connected 
with the capacity functional of the set B = {o, h}, which is 



Tu^{{o,h}) = 1 - exp |-A J^^ i^i[l\{o,7rL^z)m -K{Z)) 0(dZ)| . 

Let 7^ denote the covariogram of a measurable set A C L{Z)-^ defined by 

7A(x)=i^£?(yln(A-x)) 

for X G L{Z)^. 



(5) 
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Lemma 2. For h G M we have 



Cu^ih) = 1 - 2exp|-A J ^A{Z)e{dZ)^ 

+ exp I -2A J^^ A{Z) e{dZ) + A ^k[Z) {t^l{Z) (h)) e{dZ) | . 



(6) 



Proof. Consider the term {o,ttl(z){^)] © -K{Z) = -K{Z) U (7ri(z)(h) - K{Z)). Its 
volume is equal to 

v^,[l\-K{Z)®{o,i:Liz){m 
= 2A{Z) - v^,[1\k{Z) n {K{Z) - ^LiZ){m 



= 2A{Z)--lK{Z)i.'^L{Z){^)- 

Using equations Q and (|5|, the covariance C[7g(h) rewrites 
Cu-^{h) = 2p-Tu^{{oM) 

= l-2exp|-A j ^ A{Z)e{AZ)^ 



HZ), 

d-k ^ 



-K{Z)®{o,7TLiz){h)})e{dZ) 



+ exp|-A/^^^. 
1 - 2 exp I -A J ^ A{Z) e{dZ) | 

+ exp |-2A^^ A[Z) e{dZ) + X J^^ lK(Z){^Liz){h)) 9{dZ)\ 



□ 



Example. In the following, we give an example of a cylinder process in two dimensions 
with cylinders of constant thickness 2a where the integrals in ^ can be calculated 
explicitly. 

Let I be an arbitrary line through the origin, ip the angle between the x-axis and 
and h = {r,ip) a vector in polar coordinates. We use the notation Ba{o) x ip with 
If G [0, vr) for a cylinder with radius a and direction space I. Since |7r;(h)| = r| cos{^p—^p)\, 
formula ^ rewrites 

Cu^ih) = 1 - 2e-2^'^ + 6-4^'^+^^, 

where 

''^{2a-\7Ti{h)\)l{\7ri{h)\<2a}e{Ba{o)xdip) 

(2a - r| cos((/9 - V')!) 0{Ba{o) x d(p). 



2a 
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In the isotropic case {6{Ba{o) x dip) = dp>/n) we can choose ip arbitrarily, for example 
tp = 7r/2. This yields 

f ^ 

I = (2a — rsimp) — . 

J<pe[0,Tc]:smip<^ TT 

In case r < 2a this simplifies to 



dip 2 
I = 2a — r / sinp — = 2a — r— . 

Jo TT TT 



And for r > 2a we get 



On I /■arcsin — i-k \ 

I=-( 'd^+ ,dJ + 

J, / /-arcsin ^ rir ^ 

H / {—sirLip)dp+ / (— sin</5)d(^ 

TT yJo A-arcsin^ y 

4a /2a\ 2r / / /2a\\ 
= — arcsm I — I H (cos ( arcsm I — 11 — , 

4a /2a\ 2r / r /2a V 
= 2a arccos — 1 — Wl— — 



which gives us the final formula 

■ I _ 2e-2Xa ^ ^-2\a-^ ^ if ^ < 2a, 

Cf;,(h)=<(l-2e-''" + «^p{-2Aa- 

(^4aarccos (^) +2r (^1 - ^1 - |, if r > 2a. 

The first derivative of Cu^ (h) will be needed later for the calculation of the intensity 
iSs of the surface area measure of U-.- 

Proposition 1. Suppose that H is a simple stationary Poisson cylinder process with 
shape distribution 9 and J^o S{K{Z))9{dZ) < oo. Then the derivative of the covariance 

k 

function in direction h at the origin is given by 
C'u^{o,h) 

= Aexp |-A^^ A{Z) 0{dZ)\j^^ 7x(z)(o, vri(^)(h)) [h,L{Z)] 9{dZ), 

\ is J k 

where 7^(0,77) denotes the derivative of o-t the origin in direction r), and [^,77] is the 
volume of the parallelepiped spanned over the orthonormal bases of the linear subspaces 
of ^ and Tj. 
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Proof. To simplify the notation, we shall also write [x, ry] for t]] if ^ is the line spanned 
by X. By Q, we have 

Cusio) = p = 1 - exp |-A J ^ A{Z) e(dZ)| 



and thus 



where 



Cc/Jh) - Cu^{o) = exp |-A^^ A{Z) &(dZ)| (e^ - 1), 

J = A / j7x(z)(vrL(z)(h)) - 0(dZ). 

We observe that ^(Z) - v^[^^{K{Z) n (K(Z) - til^z){^))) is equal to zero if 
■^L(Z)(h) = o, and is less than or equal to |7r2,(^)(h)|S'(i('(Z)), otherwise. 
This yields 

|J|<A|h|/ S{K{Z))e{dZ) = 0{\\i\), h^o. 
Thus we obtain e'^ — 1 = J + o(J) = J + o(|h|) for h — o, and 



C(7^(o, h) = lim 



h-5>o Ihl 



exp{-A/^^^A(Z)«(dZ)j(lim^ 



So we need to investigate the behavior of J/|h| as h — > o. By the dominated convergence 
theorem, we get 

lim ^ = A / li,„ ^.^ir(A-(^)nWZ)-.„„(h)))-^(^)A 
h^o|h| p^kL,Z,(h)| ; 

^^£?(A-(^) n mz) - 1)) - ^(z)A ^^^^^ ,,znsm 

Jz?t^o\ \t\ / 



A^^7x(Z)(o,t) [h,L(Zp(dZ), 



where |7r/,(2)(h)|/|h| = | cos Z(h, L(Z)-'-)|, t = 7ri(^)(h), and Z{h,L{Z)-^) is the angle 
between vector h and plane L{Z)-^. □ 

3.3 Contact distribution function 

Let B be an arbitrary compact set with o €z B (called the structuring element), and 



let r > 0. The contact distribution function (cf. (Stoyan et al 1995 p. 71)) Hsir) 
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P(C/h n rB 7^ I o ^ of the union set of the stationary Poisson cyhnder process H 
with structuring element B and volume fraction p £ (0, 1) can be calculated as follows: 

^ ^ _ P{U^nrB = ill) ^ ^ _ 1-Tu-JrB) 

1 — p 1 — p 

_ ^ exp{-Xj^.^ui[l\-K{Z)®7rLiz){rB))e{dZ)} 
exp{-XJ^.^A{Z)0{dZ)} 

= 1 -exp|-A J^^ [i^i[l\-K{Z)®7:Liz){rB))-A{Z)] 0(dZ)| . 

Further simplification of this formula is possible in some special cases. 

Consider the contact distribution function Hb with B being a line segment between 
the origin and a unit vector rj. In this special case the contact distribution function is 
called linear. With a slight abuse of notation we shall use a vector to represent the line 
segment between the origin and the endpoint of the vector. It will be clear from the 
context whether the vector or the line segment is meant. 



Lemma 3. If the probability kernel (cf. is concentrated on convex bodies and 

isotropic in the first argument for all ^ E G{k, d) then for a unit vector r] the linear 
contact distribution function of U-= is given by 

Hr,{r) = 1 - e-^''Co{v) (8) 

with 



JG{k,d) JKn^^ 

Cd^k = ' '^^'^ /C n denotes the family of all convex bodies in ^-^ . 

Proof. Qlt follows from that ^ holds iff 

rCo{v) = l[v^d[l\-K{Z)®nL^z){rr,))-A{Z)\ 9{dZ). 



Using the notation introduced in (Schneider 1993, p. 275-279) for mixed volumes (here 
all mixed volumes and surface measures are w.r.t. L{Z)-^) we calculate 

v^,[l\-K{Z)®7:L^z){rv))-A[Z) 

= {d- k)V{^L(z){r^),K{Z), . . . , K{Z)) 

= ^ / \{u,T^L(Z){v))\ Sd-k-i{K{Z),du), 

where (•, •) denotes the scalar product, and Sd-k-i{K{Z), •) is the surface area measure 
of K{Z) in L{Z)-^. 



^The idea of this proof goes back to an anonymous referee. 
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Thus, 



Co{r]) = - / \{u,T:L{z){ri))\Sd-k-i{K{Z),du)e{dZ) 



K^/,7rg(7?))|5d„fe_i(K,dn)/3(dK,e)a(de). 

^ JG{k,d) J/cnc-L Js''-^n(^ 

Because of the rotation invariance of 0, the value of the integral does not change 
if we replace K with iDK for an arbitrary rotation in Furthermore, we get the 
following equation since the the surface area measure is invariant w.r.t. rotations when 
they are applied to both arguments. 

/ \{u,7T^{rj))\Sd-k-im,du)= [ \{^u,7T^{rj))\Sd-k-i{K,du). 

Thus, integration over the group rot(^-'") of rotations in ^-^ equipped with the Haar 
probability measure leads to 

\{u,TT^{r]))\ Sd^k-i{K, du) 



\{u,7T^{ri))\Sd-k-i{K,du) d'd 



\{^u,7T^{r]))\d^Sd-k-iiK,du) 
2cd,k SiK)[C,v], 



where Cd,k is the constant from the claim, and we used (Spodarev, 2002 Corollary 5.2) 
for the last equality. 
This leads to 



□ 



Co{v)=Cd,k / S{K)P{dK,0[LvHdO- 

JG{k,d) Jicni^ 



Now let the structuring element B be the ball Bi{o). In this case the contact distri- 
bution function is called spherical. It is obvious that ■Ki(^2:){Br{o)) is a ball of radius r 
in the {d — A:)-dimensional subspace L{Z)-^. If K{Z) is almost surely convex then the 
use of the classical Steiner formula leads to 



d-k 



Eu^^i\-K{Z) © nLiz){Br{o))) = EA{Z) + ^ K,EVtlmZ))r\ 



1=1 

which yields 

d-k 



HB,io){r) = 1 - exp -A ^ f ^ Vtl,{K{Z)) e{dZ) . 

{ i=i -'^k ) 
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Example. In what follows, the case of dimensions two and three is considered in detail. 
It is assumed that the conditions of Lemma [3] hold. 

• For d = 2,k = 1 Lemma |3] yields 

JG(1,2) Jjcn^^ Jg{i,2) 

Hence, it holds Hn{r) = 1 — exp |— 2A r Jq^^ 2) [^^ "H] Q^(dC)|) and so Hjj{r) does not 
depend on K{Z). 

And for the structuring element being B = i?i(o) one gets 



HB,{o){r) = 1 -exp|-2Ar^^yo'(^(^))^(d^)| = 1 



e-2^^ 



Interestingly the result does not depend on the distribution of the cross section. 
For d = 3,k = 1 we get 

Co(r?) = - / / SiK)PidK,Cm^HdO 
TT JG(i,3) JKn^^ 

which yields 

H^{r) = l-exp\-— I I S{K)P{dK,m,r,]a{di) \ . 
vr JGii.-i) JlCnS^ 



For K{Z) = Ba{o) we have 



vr 

Thus, 



G{k,d) 



And if the structuring element is the unit ball {B = Bi(o)) then 

HB,(o){r) = 1 - exp |-A (^2r J^^ V^\K{Z)) e{dZ) + J^^ K2 e{dZ) 

= 1 -exp|-A l^r^ S{K{Z))e{dZ) +r\^'^ , 

where S{K{Z)) is the perimeter of K{Z). 

If additionally K{Z) is a ball of constant radius a then 

TT (^\ 1 27raAT-— TiAr^ 
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4 Specific surface area 



In the recent paper Hoffmann ( 2009 1 , the specific intrinsic volumes of a rather general 
non-stationary cylinder process are given. In the stationary anisotropic case, some of 
these formulae can be simplified. In this section, we give an alternative proof for the 
specific surface area of the union set Us of a simple stationary anisotropic Poisson 
cylinder process H leading to a simpler formula than that of Hoffmann (20091 which can 
be immediately used in applications. 

The specific surface area 5= is defined as the mean surface area of [/= per unit volume. 
More formally, consider the measure Sus{B) = EW^-^idUs n B) for all Borel sets 
B C W^. We assume that this measure is locally finite, i.e. Su^{B) < oo for all compact 
B. Sufficient conditions for this can be found in Lemma |4j Due to the stationarity of 
H, the measure Su^ is translation invariant. By Haar's lemma, there exists a constant 
Ss ^ such that Suv,{B) = S'^Ud^B) for all Borel sets B, cf. Ambartzumian (19901. 
The factor 5= is called the specific surface area of [/=• 

Lemma 4. The specific surface area 5= of the union set {/= of a stationary anisotropic 
cylinder process E is finite if S{K{Z))6{dZ) < 00. 

Proof. Let B := Bi{o) be the unit ball about the origin. Then we calculate using the 
abbreviation Lq = L{Zo) and Campbell's theorem 

Su^(B) = ETi'^-^idUs nB)<EY^ n'^'HdZ D B) = f n'^-\dZ D B)A{dZ) 

/ / n''-\{dZo + x) n i^tki^^^) ^(d^o) 

[ [ [ lB{y)n'-Hdy)i.ji^,{d^)9{dZo) 

f f f lB(y + x)?^'^-i(dy)z.,^_!,(dx)0(dZ,) 
Jz° Jl^ JdZo 

(vru(y))l.,jB)(vrL„(y) + x)i.f_!,(dx) T^-^-^Cdy) 0(dZo) 



= A 
= A 
= A 
< A 
= A 



Z-JdZoJL^ 



Z° JdZo 



B)i7rL^iy))u^l,{7rLom'H''-\dy)eidZo 



d-i, 



X'^tki^LoiB)) I n^'-HdZoniTTL^iB) X Lj))e{dZo) 
I i^ii,i7TL.iBm'-'-HdKiZo))eidZo) 



S{K{Zo))B{dZ,). 



This yields 5= = Sus{B)/ud{B) < 00. 



□ 
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The following results hold for any random closed set X with realizations almost surely 
from the extended convex ring S which is defined as the family of sets B with BnW £ TZ 
for any convex compact observation window W. 

Lemma 5. Let X £ S be an arbitrary stationary random closed set with finite specific 
surface area. Then the specific surface area of X is given by 



Sx 



dKd 



l^d-l JG{l,d) 



A(e)de, 



(9) 



where d^ is the Haar probability measure on G{1, d), = \E'^q{X n ^, -Bi(o) n ^) is 
the intensity of the number of connected components of X on a line ^ G G{\,d). 



Proof. By Crofton's formula for polyconvex sets (cf. (Schneider and Weil 2008 Th. 6.4.3)) 
and Fubini's theorem, we have 



Sx 



1 

2T{ 



En'^-^{dxr\Bi{o)) 



Kd 



-E^a-i{X,Bi{o)) 



d+l- 
2 ' 



KdT{d/2) 
dKd 



e[ f ci>o(xn(c + x),i3i(o)n(e + x))i.«_,(dx)de 

Jg(14) M-^ 



1 



i^d-1 Jg(i4) 2 



/G(i,d) Ji 
i?a>o(xnc,i?i(o)nOde 



The following result generalizes the well-known formula 

dud 



Sx 



K-d-l 



-C'xiO) 



□ 



(10) 



(Stoyan et al 1995 p. 204), for stationary, isotropic, and a.s. regular random closed sets 
X € S to the anisotropic case. A closed set is called regular if it coincides with the 
closure of its interior. Note that, since in the isotropic case Cx{h), /i E M'^ depends only 
on the length of h, and not on h itself, in this formula Cx is a function of a real variable, 
namely the length of h. For the particular case of stationary anisotropic random sets in 



formula (11) can also be found (without a rigorous proof) in Berryman (1987). 



Theorem 1. Let X be an a.s. regular stationary random closed set with realizations 
from S and finite specific surface area. If Cx(h) is its covariance function then the 
specific surface area of X is given by the formula 



dKd 



Kd-l JG{l,d) 



C'x{o,v^)di, 



(11) 



where C^(h, v) is the derivative of Cx(h) at h in direction of unit vector v, and is 
a direction unit vector of a line ^ G G{\,d). 
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Proof. For a stationary random closed set C/ C M from the extended convex ring denote 
hy —U the set reflected at the origin. Define a random variable V which is uniformly 
distributed on {—1,1} and independent of U. The random closed set UV is obviously 
isotropic, and thus formula (10) yields Suv = ~2C^y(0). Since Su = Suv and C'uifi) = 
C[;y(0), this means that Su = -2C[^(0). 

Applying this to ?7 = X n e, e G G{l,d), we get A(e) = \Sxn^ = -Cxnc(O) = 
— C^(o,rg). Lemma [5] completes the proof. □ 

If X is an a.s. regular two-dimensional stationary random closed set with realizations 
in 5, formula (111 simplifies to 

Jo vr JO 

The following result is a direct corollary of Proposition [T] Theorem [T| and Fubini's 
theorem. 

Corollary 1. Let E be a stationary Poisson cylinder process with intensity X, shape 
distribution 9 and cylinders with regular cross-section K{Z) S IZ for 6-almost all Z ^ 
and finite specific surface area. Then, the specific surface area of Us is given by the 
formula 

Ss=-X^[ I 7^(z)(o,vr^(z)(r5)) K,M^)]de0(dZ)x 



X exp I -A j A{Z) e{dZ) 



Example. Assume that K{Z) is convex and regular for 0-almost all Z G Z^. 
• For arbitrary d, and k = d — 1 \i holds for d^-a.e. line ^ G G(l, d) that 

7^{z)(o,vri(z)(r5)) = -1 

and 

(d+l)Kd+iKi {d+l)Kd+l 



j [tL{Z)]d^= f [e,L{Z)]d^ 

JGd.d) JG(d-l.d) 



G(l,d) JGid-lA) dKd2K2 dUdlT 



see ( |Spodarev| [2002| Corollary 5.2). 
This yields 

{d+ l)Kd+l J f L{Z) 



= ;, (^+ 1)^^+1 exp -A / u^^'\K{Z))e{dZ) 
= 2Aexp I -A J v[^^\k{Z)) e{dZ) I . 
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For d = 3, k = 1, K = Ba{o) it can be calculated that 7^(2)(o, 7ri(2)(^)) 

/g{i 3)['^' ^i^)] dC = 1/2 (see also (Stoyan et al 
Corollary 5.2)), and thus we have 



-vra, 



1995 



p. 298), or (Spodarev 



2002 



4X-7rae 
2 



27raAe" 



-Xna-' 



which coincides with the case of isotropic cylinders, compare (Ohser and Miicklich 
2000] p. 64). 



5 Optimization Example 

In this section we show how the formulae from Sections [3] and |4] can be applied to solve 
an optimization problem for cylinder processes. 

The following problem originates from the fuel cell research. The gas diffusion layer 
of a polymer electrolyte membrane fuel cell is a porous material made of polymer fibers 
(see Figure [T]) which can be modeled well by an anisotropic Poisson process of cylinders 
in M^. In a gas diffusion layer, the volume fraction of the polymer material lies between 
70 and 80 percent, and the directional distribution of fibers is concentrated on a small 
neighborhood of a great circle of a unit sphere 5^, i.e. all fibers are almost horizontal. 
In order to optimize the water and gas transport properties, it is desirable to have a 
relatively small variation of the size of pores in the medium, where we define a pore at 
a point X in the complement of Us as the maximal ball with center in x which does not 
hit Ue. 

We investigate the following mathematical simplification of this problem, which can 
be solved analytically in some particular cases. 

For a fixed intensity A of the Poisson cylinder process H, find a shape distribution of 
cylinders 6 which maximizes the volume fraction p of [/= provided that the variance of 
the typical pore radius H is small. In other words, solve the optimization problem 

ip ^ max,, 
[VarH < e, 

where H is a random variable with distribution function H^^^^^-j^r). 

As it will be clear later, the condition on the directional distribution a of fibers that 
all fibers are almost horizontal can be neglected since the directional component of the 
shape distribution 9 has no influence on the solution. 

To simplify the notation, let Cs = J^o S{K{Z))9{dZ) and <I>(x) be the distribution 
function of a standard normally distributed random variable. 

First we take a look at the moments of the pore radius H (assuming that r > 0) 
remembering that /^^^(^^(r) = 1 — e~-*'('''^''+'~ (as shown in an example in Section 



3.3) 
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and thus the density of H is equal to ^Hbj^(o){^) = K^s + 27rr)e Hrcs+r'^Tr) _ j^. j-^Q^jg 
EH = rX{cs + 27rr) exp ^-vrA (^r + ^ exp j dr 



oo 



47r i ^ 



Furthermore it can be calculated that 



EH'^ = exp 



Defining Cg = exp (4^) and c$ = ^1 — <I> y^sy^jj, this leads to 

'''' "^''''^ WT" — -"' 

multiplication with vrA yields the equivalent condition 

1 — \/AceC$Cs — 7rCgc|, < evrA, 

which holds if and only if 

C\/Acc\ Ac^ 
cec^ + ^j -^ + (eA-l/vr)>0. 

This is always fulfilled if e > and — (eA — l/vr) < or, equivalently, < 

In the following we always assume that £ > and replace the condition Vari? < e 
by a stronger sufficient condition 



J^^S{K{Z))9{dZ) <27r^e-^^. 



(13) 



Hence, (12) is reduced to the optimization problem 

/^o A{Z)d{dZ) — maxg. 



J^oSiKiZmdZ) <27r^e-^,. 



(14) 
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The solution of the optimization problem ( 14 ) yields cylinders with 6-a.s. circular base 



Notice that this solution does not depend on the directional distribution component a 
of 9. Indeed, cylinders Z can be replaced by cylinders Z' which have the same direction 
space and surface area {S{K{Z)) = S{K{Z'))) but are circular. Then the isoperimetric 
inequality yields A{Z') > A{Z). Thus, it holds that 

S[K{Z))9{dZ) = [ S{K{Z'))9{dZ) 
f A{Z)9{dZ) < f A{Z')9{dZ), 

Jz° J2° 

which means that the circular version is at least not worse than the original version. 
Thus, we assume that the cylinders are 9-a.s. circular and denote the radius of a 



and 



cylinder Z by R{Z). It follows from condition (13) that 



/ ^ S{K{Z))9{dZ) =27: I R{Z)9{dZ) < 2n^e - ^, 



i.e. the new condition is that the expectation of the radius of a typical cylinder is less or 
equal than ^ 



ttA- 



Furthermore, it follows from (14) that maximizing p is equivalent to maximizing 
j^oR{Zf9{dZ). 

The above calculation shows that the volume fraction of 70% — 80% in the optimized 
gas diffusion layer of a fuel cell can be achieved best by taking fibers with circular cross 
sections, relatively small mean radius and high variance of this radius. 

Figure |2] shows that cross sections of fibers of gas diffusion layers are almost circular. 
There are also gas diffusion layers with a little variance in the fiber radii, although they 
are mostly nearly constant. Anyhow the variance of the fiber radii is of course limited, 
since it is impossible to produce fibers with an arbitrarily large radius. 



We have to remark that from a practical point of view the optimization problem ( 12 ) 
is not well posed. For the construction of gas diffusion layers, mainly the intensity of 
the fibers A can be varied. Hence a practically relevant optimization should involve 
maximizing the volume fraction p with respect to A as well. Since the latter problem is 
much more involved than the one discussed here, it would go beyond the scope of this 
paper. 
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' 300 /im ' 

Figure 2: Microscopic picture of the gas diffusion layer 
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